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Abstract. Quarkonia are the central objects to explore the non-perturbative nature of non-abelian 
gauge theories. We describe the confinement-deconfinement phases for heavy quarkonia in a hot QCD 
medium and thereby the statistical nature of the inter-quark forces. In the sense of one-loop quantum 
' effects, we propose that the "quantum" nature of quark matters follows directly from the thermody- 

namic consideration of Richardson potential. Thereby we gain an understanding of the formation of hot 
and dense states of quark gluon plasma matter in heavy ion collisions and the early universe. In the case 
of the non-abelian theory, the consideration of the Sudhakov form factor turns out to be an efflcient tool 
for soft gluons. In the limit of the Block-Nordsieck resummation, the strong coupling obtained from the 
Sudhakov form factor yields the statistical nature of hadronic bound states, e.g. kaons and Ds particles. 
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1. Introduction 

^ . In this work, we study the geometric nature of the quark matter formation. Specifically, we shall 
^ I illustrate that the components of the vacuum fluctuations define a set of local pair correlations against 
the vacuum parameters, e.g. charge, mass and angular momentum. Our consideration follows from the 
notion of the thermodynamic geometry [TH3]. Importantly, this framework provides a mathematical 
platform to exactly understand the nature of the pair local correlations and underlying geometric 
structures pertaining to the global phase transitions in quarkonium systems. This perspective yields 
a well-known understanding for the phase structures of mixtures of gases, black holes in string 
theory [3H7j and in other diverse contexts, as well. 

The main purpose of the present investigation is to determine the thermodynamic properties of 
the quarkonium configurations, in general. Quantum chromodynamics (QCD), as the theory of strong 
interactions, celebrates physics [SHlOj at both the high and low temperature domains. Thereby, our 
consideration plays a crucial role in understanding the phases and stability of the matter formation. 
In the soft gluon limit, viz., for a small transverse momentum k±, the behavior of the abelian theory 
follows directly from the Poisson distribution and its numerical counterparts. On the other hand. 



the statistical nature of the non-abehan soft gluons is understood in terms of the Sudhakov form 
factor |12j . Let us recall that the QCD coupling as{k±) never lies near the limit k± — t- 0, and so the 
QCD effects are limited for the bound state thus formed after the resummation. As mentioned in the 
Ref. |15j . an integration over k± requires the Sudhakov form factor 

h{b) = (1) 

where k± € (0,'n^p/2) is due to a physical reason, i.e. that on average each soft gluon can take 
as much as half of the initial center of mass energy. Towards the determination of the index p, an 
interesting argument follows from the work of Polyakov [11] which we shall explore further from the 
perspective of thermodynamic geometry in the subsequent consideration. Before doing so, let us 
consider the joint effects of the (i) confinement and (ii) rotation, and thus make a platform to describe 
the thermodynamic geometry of rotating quarkonia. Such a simplest configuration is described by 
Regge trajectories with the leading order effective potential 

V{r,J) = :^^^^-p^ + Cr^P-\ (2) 

In this case, the effective theory |13lll4j is inspired from the limiting QCD strong coupling 

QCn 

where b := (33 — 2A^j)/127r. In the sense of one-loop quantum effects, we propose that the "quantum" 
nature of quark matters follows directly from the thermodynamic intrinsic geometry of the Richardson 
type potentials. In the momentum space, the net effective potential offers the right quarkonium bound 
states, after taking account of the one-loop exchange terms. 

Let us now focus our attention on the thermodynamic geometry of quarkonium bound states with 
finitely many parameters of the effective field theory. We consider here a framework of the intrinsic 
Riemannian geometry whose covariant metric tensor is defined as the Hessian matrix of the QCD 
coupling, with respect to a finite number of arbitrary parameters carried by the soft gluons and 
quarks. Such a consideration yields the space spanned by n parameters of the strong QCD coupling 
as, which, in the present treatment, exhibits a n-dimensional intrinsic Riemannian manifold M„. As 
per the definition of the thermodynamic geometry [IH7] , the components of the covariant metric tensor 
are given by 

■" dx^dx' ' ^ ^ 

where the vector x E M„. In the strongly coupled quarkonium effective configuration, there are only 
a few physical parameters, what makes the analysis fairly simple. As mentioned in the Ref. [15], the 
variables of the interest of the above quarkonium configurations are the momentum scale parameter, 
:= q, the mass M, and the angular momentum J, if any. 



2. Massless Quarkonia 

Let us first examine the massless non-rotating quarkonia and study the thermodynamic stability 
properties as emphasized earlier, and subsequently include the rotation. Considering Eqn.(l3]), one 
obtains the following expression 

61n(l +p{q/L)P) 

for the strong QCD coupling, where L := ^qqd- '^'^ compute the thermodynamic metric tensor in 
the space of the parameters {q,p}, we employ the Eqn.Q, which leads to the following expression for 
the components of the metric tensor 



In order to simplify the subsequent notations, let us define the logarithmic factor as 



lip):=Hl+p[g/Lr). (7) 

Thus, the factors in the numerator of the local pair correlations can be expressed as 

2p2(<j/L)2P + l{p){{q/Lr - Piq/LY + {q/hf^p), 

2{q/Lf^p + lip){2{q/LfPp2 _ g^^/^^p _ 2{q/LfPp - (g/L)fpln(9/L)), 

2(g/L)2P(p + 2p2 ln(g/L) + \n{q/Lf) - l{p){{q/LY + [q/lf^'p 

+A{q/LYp\n{q/L) + 2{q/Lf^p^ \n(q/L) + (q/LYp^ Hq/hf). (8) 
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In the case of the massless non-rotating quarkonia, we find that the factors {r^| i = 1, 2, 3} take the 
uniform value exp (2/(p)) in all denominators of the local pair correlation functions. In a given 
QCD phase, this happens when the parameters {q,p} are confined in the domain 

{V := {q,p) e A/2 1 nf^ > 0, ng, > 0}. (9) 

Over the above domain of {q,p}, the massless non-rotating quarkonia is well-behaved and locally 
stable. From the definition of the thermodynamic geometry, we find further for the generic value 
of the parameters that the Gaussian fiuctuations form a stable set of correlations over {q,p}, if the 
determinant of the metric tensor 

remains a positive function on the intrinsic g'p-surface {M2{R), g). Explicitly, we obtain that the 
numerator of the determinant of the metric tensor can be expressed as 

nf[q,p): = 2{q/ Lf^ {p + Zp'^ + 2p^ln{q/ L) + 2p^ln{q/ L) + p^ln{q/ Lf) 

-l(p){2(q/LfP +p{q/Lf^ + 2p\q/Lfnn{q/L) + Ap{q/Lfnn{q/L) 

+pHq/L?''ln{q/LY + 2p' {q / LfHniq / L) + ApHq/L?" + Ml/L?")- (H) 

In this case, we find that the scalar curvature reduces to the following specific form 

We find that the factors of the numerator of the scalar curvature take a set of interesting expressions 
|15] . In the case when the Ricci scalar curvature R{q,p) vanishes, the underlying quarkonium system is 
found to be in equilibrium. Such a state of the configuration can arise with {nW'S = 0, i = 0, 1, 2, 3}, if 
the other factors of the scalar curvature remain non-zero. In the other case, when the scalar curvature 
R{q, p) diverges, the configuration goes over a transition. Such an extreme behavior of the quarkonia 
is expected to happen, when either the index p or the numerator of the determinant of the metric 
tensor vanish. We further observe that the stability of massless non-rotating quarkonia exists in 
certain bands. A closer view shows that the quarkonia are decaying and interacting particles in the 
Coulombic limit. For the regime of the rising potential, the respective determinant of the metric tensor 
and scalar curvature shows that the limiting rising potential quarkonia are stable and non-interacting 
particles. In the Regge trajectory model [15], we find interestingly that all possible local and global 
thermodynamic stability behavior of the three parameter {q,p,J} massless quarkonia remains the 
same up to the sign of bi, as if there were no effects of the rotation in the underlying configuration. 
The fact that the efficiency of the rotation induces a mass to the quarkonium is analyzed by considering 
the Bloch-Nordsieck resummation of the angular phases. 



3. Massive Quarkonia 

Let us firstly illustrate the cases for the two parameter configurations with either {q, J} or {q, m} 
fluctuating and then systematically extend the consideration for the generic three parameter quarkonia. 



3.1. Two Parameter Quarkonia 

For the of massive quarkonia, the resummed strong QCD couphng takes the following form 



' b ln{l+p{q/L)P) ^^-^rj^'' ^ > 

where J(i^, x) is the Bessel function of the first kind of the order i/. In order to describe fluctuations 
in the QJ-plane, let the logarithmic factor concerning the Bloch-Nordsieck rotation be defined as 

/(,,J) :=ln( ^ + ^ ). (14) 

From the Eqn.(j3]), we find that the components of the metric tensor are 

p{n^°,^-' + n[\^'l{p) + nf^^-'l(pf) pjnf^^' + n\\^' l{p)) p(2J - q){l - J(0,a^)) 

4M{p)3exp(2Z{p))gV2(j_q)3/2 • ^'-'P 2M{p)2 cxp (/{p))g3/2 v^(J - g)3/2 ' 2W{p) j3/2( J _ g)3/2 ' ^ 

Without any approximation, the factors in the numerator of the pure qq and (7 J-components can be 
expressed as linear combinations of the integer powers of the scaling (g/L)"^, where n ^ Z. As a 
result, we see that the geometric nature of the parametric pair correlation functions turns out to be 
remarkably interesting, viz., the domains of the local stability of the fluctuating quarkonia may be 
easily described in terms of the parameter q and J. Under the Gaussian fluctuations of {q, J}, the 
local stability of the system requires that (i) qq- fluctuations satisfy the constraint 

nfl' + n'-^^h{p)+nfyi{pf<0 (16) 

and (ii) J J- fluctuations be constrained to the following limiting values of the Bessel function 

J{0,a^) < 1, q>2J, 

> 1, q< 2J. (17) 

In this case, we flnd that the determinant of the metric tensor reduces to the following expression 

|, ||_ p\4'^-'+n^^^'l(p) + nf^'l(pr) 

8b2i(p)4exp(2«(p))g5/2j3/2(j_^)3/2' 

where 4^^"^ = 4p\q/L)P{n^g^'>-^ + n'g'^^-^p{q/L)P) and = °^^ + 2nf ^^■^p(g/L)P + nf ^^^(<?/L)2V- 
Furthermore, it turns out that all factors appearing in the numerator of the determinant of the metric 

tensor, e.g., {ng*^^"^, n^^^^"', n^^'^^'^, n^^*^^"^, n^^"'^^"', ng^"^^"^}, can be presented as linear combinations of the 
scaling (q'/L)"^, where n £ Z. Combining the effects of all fluctuations of the {q,J}, we observe 
that the quarkonia are stable for q := G (1)4). In general, the global stability requires that the 
determinant of the metric tensor must be positive deflnite, which in the present case transform as 

n^g -.= 4'^^ + 4''>-'lip) + n'-P''lipf < 0. (19) 

It turns out that the thermodynamic curvature may be written as the series of the charmonium 
logarithmic factor l{p) and Bloch-Nordsieck logarithmic factor f{q,J) of rotation as the coefficient of 
the expansion. Systematically, the exact expression for the scalar curvature takes the form 

where the Bn in the numerator of the scalar curvature are polynomials in p, whose coefficients are the 
functions of the Bloch-Nordsieck logarithmic factor f{q, J). In sequel, we find that the quantitative 
properties of the scalar curvature and the Riemann curvature tensor remain similar |15j . when we 
consider the variable Q as the transverse momentum with the understanding that k = k± and the 
mass as the other variable of the arbitrary two parameter real quarkonium system. The present 
consideration explicates the regions of the thermodynamic (in) stability for the massive quarkonia and 
the (un)stable phases of the one loop QCD. As the non-linear effects become stronger and stronger, 
it turns out that the thermodynamic instability and correlations grow further. This motivates us to 
extend our analysis to the case of more general quarkonium configurations. 



3.2. Three Parameter Quarkonia 

Finally, let us analyze the general massive rotating quarkonia, when all parameters, viz., the scale q, 
index p and angular momentum J of the theory are allowed to fluctuate. In the framework of the 
Bloch-Nordsieck resummation, the strong QCD coupling takes the following form 

' b ln{l+p{q/L)P) ^^+^rjzr^> ^ > 

After some simplification, we find that the components of the metric tensor are 

4W(p)3exp(2/(p))g5/2(j_q)3/2 ' ^IP 2W{p)3 exp (2/(p))q3/2 ( J - q)l/2 ' 

p(n^°3'° + nW^;(p)) _ /(-I + J(0, a7g))(4°'^ + ng^'^lip)) 

^"■^ 2W(p)2exp(i(p))g3/2(J-q)3/2ji/2' W(p)3 exp (2/(p)) 

(-1 + J(0,av/g))(4°3^'^ +;(p)nW°) _ p(-l + J(Q,a^{qm2J - q) 

^"■^ bl{p)^exp{l{p)){J -qy/'^J^/'^ ' 2b/(p)(J-g)3/2j3/2 ' > 

where the coefficients {nj^^^*^, n^^^*^, n^j^*^, n^g^*^} appearing in the components of the metric tensor 
are shown to factorize as before. In the powers of p, the exact factorizations are given as follows: 
(i) the gq-component 

= 4p2{g/L)-n(f)«+4p3(,/L)2Pn(f'«, 

= n^f)«+2p(,/Lrn(f)«+p^(5/L)^^ (23) 

(ii) the gp-component 

= n(f« + 2p(,/L)^nf,^)«+p^(,/L)2^nr° (24) 

Interestingly, the factors appearing in the numerator of the qq, qp and g J-components can be cascaded 
as a linear combination in the integer powers of the scaling (q/L)^^, where n € Z. Similarly, the 
corresponding factors of numerator of the pp-components take the following expressions 

4°)'^ = -2p{q/LfP{l + 2pln{q/L)+p^ln{q/Lf), 

4V'' = (<?/L)f(2 + 4pln(q/L)+p2l„(g/L)2)+p(q/i)2p(i + 2pln(q/L)). (25) 

Finally, the factors in the numerator of the pj-components are 

^23"" = -P(q/LY -p'Mq/L){q/LY, 

= 1+P{q/LY. (26) 

The local stability of the configuration requires a set of constraints on the domains of the parameters. 
Specifically, for the same sign of {p, 6}, we find that the local stability of the fluctuating quarkonia 
enforces the following simultaneous requirements 

(i) the qq- fluctuations satisfy 

n(r + n««^(p)+n(r/(p)2<0, (27) 

(ii) the J J- fluctuations remain within the limiting values of the Bessel function 

J(0,av^) > 1, 2J <q, 

< 1, 2J > q, (28) 



(iii) the pp- fluctuations satisfy 



nf.^° + tS^.^°1{p) > 0, J(0,av^)>l, 

< 0, J(0, a^) > 1 (29) 



for the same sign of {f,b}. The silent feature of the JJ-fluctuations is that we find the two distinct 
local behaviors for J > and J < 0. Subsequently, we find that the gp-surface is stable, if there exits 
a positive surface minor 



^•^ ' 462Z(p)5exp(3Z(p))(j5/2(J_g)3/2 • ^ ' 

It is worth mentioning that {n^g^^ , n^p'^ , n^^'^ , n^g^'^ } have the following factorizations 

= (g/L)3f„(?3)G, 
n(?)« = (g/L^nf )^ + (q/L)2fnf + (q/L)3fnf' °, 

In this identification, the factors appearing in the powers of {q/L)P have the following structures 
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(32) 



Further, the factors 032, 12i, 13i, 21i, 22i and 23i take the following explicit expressions 
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(33) 

Interestingly, the sub-factors of the factors 032, 120, 121, 122, 130, 131, 211, 212, 220, 221, 222, 
230, 231 and 232 can be expressed as the linear combinations of the integral powers of the scaling 
ln{q/L) and the integral powers of the zeroth order Bessel function. As per the expectation, we observe 
that the 3ij-factors that appear are relatively simpler than their foregoing counterparts. In fact, we 
find that the 30-factors can be seen as a two variable polynomial in J, q. For the general quarkonia 
with {q,p, J} fluctuating, the stability of the qp-suviace requires that the principle minor remains 
positive on (M^jq). Correspondingly, this leads to the constraint that the {q,p, J} satisfy 

n(0)« + nW°i(p)+n(?)«/{p)2+„(^)«/{p)3 < 0. (34) 

As per the generic fluctuations of the strongly coupled massive rotating quarkonia, we obtain the 
following general expressions for the determinant of the metric tensor 



P/(-l + J(0."V^)) ,JO)G ^ „(1)G„„, ^ „(2)G„^,2 ^ „(3)G 

863;(p)6cxp(3/(p))g5/2j3/2(J- g)3/2 



(4"'' + 4^^"i(p) + 4'>^lipf + n^/^^lipn (35) 



As mentioned before, we find that the {n^'''*'^, ng^"**^, n^^'*'^, n^^^'^j-terms factorize as follows 

4°)^ = (9/L)3^'(8p4rrr'' + 8p^nf)« + 8p«nf'«), 
4i)« = (,/L)2-n(^^)° + (,/L)3^'4^^)«, 

^(3)G ^ 430)G^(^/^),^„(31)G^(^/^)2p„(32)G^(^/^)3p433)G_ (3^^ 



After a direct simplification, we obtain the sub-factorizations 
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(37) 



As per our computation, we observe that the terms appearing in the various powers of l{p) have a sub- 
factorization in the index p. Specifically, we find that the factors of 12, 13, 21, 22 and 23-components 
are given by 
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As in the case of the surface minor, we observe that all the individual sub-factorizations, (e.g. 040, 
041, 050, 051, 060, 061, 120i, 121j, 130i, 130j, 2303, 2304, 2313, 2314), can be expressed as the linear 
combination over the integral powers of the zeroth order Bessel function and the scaling In [q/L)^ . In 
the case of the l{pY terms, we find that the sub-factors pertaining to 30, 31 and 32-factors have no 
dependence on In (q/L)^ and thus they are expressible as polynomial expressions in the Bessel function 
only. From this observation, we predict that the regions of the thermodynamic stability are present 
for g, J G (1,4). Globally, the stability of (Ma,^) constrains the principle minors {5m,P5 , ||5||} to 
remain positive. Specifically, for the same sign of {b,p, /}, the volume stability of the (M3, (7) imposes 
the following constraint 

3 

^ni^^°l{py > 0, J(0,av^) < 1, 
i=0 

< 0, J{0,ay^) < 1. (39) 

Importantly, it is worth mentioning that both the limiting configurations with J = q and J(0, ayfq) = 1 
are abided from the thermodynamic stability constraints. To summarize the phases of generic 
quarkonia, the exact formula for the scalar curvature may analogously be deduced, as the one we 
have offered for the fluctuations in the gJ-plane. In this sense, we find that the summation over 
l{p) naturally arises with the i?„ as the polynomials in p, whose coefficients can be expressed as the 
functions of the Bloch-Nordsieck logarithmic factor /(g, J). Based on the analysis of the present paper, 
the thermodynamically stable index of the Bloch-Nordsieck rotating massive quarkonia is constrained 
by the following set 

Vm^O ■■= {P I > 0, nf > 0, > 0}. (40) 

As the gluons become softer and softer, we find, in the limit of Bloch-Nordsieck resummation, that the 
underlying Sudhakov form factor offers all possible thermodynamically stable phases of the strongly 
coupled quarkonia. As examined for the Q J-plane, we would like to explicitly understand properties 
of the set Vm^o and the associated Ricci scalar curvature. Up to a phase of QCD, the global properties 
of the three parameter quarkonia remain the same as we have exactly indicated for the two parameter 
quarkonia. As mentioned in the Ref. |15j . we have computed the intrinsic geometric properties of the 
Bessel function of the first kind convoluted with two logarithmic functions of the respective weights 
(0, 2p). Further analysis of the geometric features of these exploitations is left for the future. 



4. Conclusion and Outlook 

We have examined tiie role of the thermodynamic intrinsic geometry for a class of quarkonium 
configurations. We have offered a geometric perspective to the confinement- deconfinement phase 
of (heavy) quarkonia in a hot QCD medium and thereby described the statistical nature of the inter- 
quark forces. Specifically, the intrinsic geometric analysis provides a set of physical indications encoded 
in the geometric quantities, e.g., the scalar curvature and possible geometrically non-trivial invariants 
offer the global correlation properties of an ensemble or subensemble of equilibrium configurations. 
In the sense of statistical mechanics, our analysis involving a Gaussian distribution of the particles 
ensures the thermodynamical properties of the underlying quarkonia in the late time limit. From the 
perspective of one-loop quantum effects, the nature of quark matter is shown to follow directly from 
the thermodynamic consideration of the Richardson potential. 

Our study of the quarkonia could further be explored towards other configurations concerning the 
non-perturbative and non-abelian nature of the gauge theories. Such a consideration provides a unified 
description, encompassing all the regimes of QCD at finite temperature, i.e. the Coulombic, the linear 
rising and the Regge rotating regimes, for both massless and massive quarkonia. Phenomenologically, 
our results can be thus be used to investigate the statistical nature of soft gluons and the associated 
phenomenon at the LHC. 
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